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O ! Abstract 



We report on surface effects on the electronic properties of interfaces in epitaxial LaAlOs/SrTiOa 
heterostructures. Our results are based on first-principles electronic structure calculations for well- 
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\^ \ relaxed multilayer configurations, terminated by an ultrathin LaAlOs surface layer. On varying the 



thickness of this layer, we find that the interface conduction states are subject to almost rigid band 
"^ ■ shifts due to a modified Fermi energy. Confirming experimental data, the electronic properties of 

heterointerfaces therefore can be tuned systematically by alterating the surface-interface distance. 
We expect that this mechanism is very general and applies to most oxide heterostructures. 

PACS numbers: 73.20.-r, 73.20.At, 73.40.Kp 
^ ' Keywords: density functional theory, surface, interface, SrTiOa, LaAlOa 
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FIG. 1: (Color online) LaAlOa/SrTiOa (001) heterostructure terminated by an (A102)^ surface, 
where the bulk is formed by alternating (Ti02)'^ and (SrO)" layers. The band structure calculations 
presented in this paper refer to surface layers comprising 3 and 4 LaAlOs unit cells. 



Heterostructures based on II-IV semiconductors have attracted erreat attention in recent 
yea. due to a la.,e ...be. of possible app„ca«o„s fo. e.ec.on.c devees Q. Special 
interest has focussed on perovskite heterostructures from transition metal oxides, because 
of both the variability of the perovskite structure and the correlated electron properties 
associated with transition metals [2]. Pulsed laser deposition techniques and molecular 
beam epitaxy nowadays allow to grow layered structures with a precision of one unit cell. It 
consequently is possible to create atomically sharp interfaces between perovskite oxides with 
strongly different electronic features, opening a huge field of new cooperative phenomena. 
Charge redistribution plays the key role for the electronic structure at the interface since 
electrons flow from one layer to another, driven by different electrochemical potentials in the 
component materials [Sj. The physical properties of such a multilayer structure hence may 
not be present in either of its constituents. For instance, a conducting quasi two-dimensional 
electron gas is found to be formed between the Mott insulator LaTiOs and the band insulator 
SrTiOg [41. 

A contact between two common band insulators with the large band gaps 5.6 eV and 
3.2 eV is realized in the LaAIOs/SrTiOs heterointerface. The structure consists of alternating 
(SrO)°, (Ti02)°, (LaO)+, and (A102)~ layers, see the schematic representation in figuredl At 
the electron-doped (Ti02)''/(LaO)^ contact a quasi two-dimensional electron gas with very 
high carrier density jsl, 61] is formed. In contrast to the mixed valence system LaTiOs/SrTiOs, 
the formation of the electron gas here is considered to result from a polarity discontinuity. 
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FIG. 2: (Color online) Band structure and partial Ti 3d and O 2p DOS (per atomic site) for bulk 
SrTiOs- The weighted bands highlight contributions of the Ti 3d t2g states. 

Because of different formal valences of the metal ions Ti^"*" and La^"*", uncompensated cliarge 
is left at the interface, giving rise to the metallic layer. Oxygen vacancies introduced in the 
SrTiOa layer during the deposition process likewise are a source of doping Tj. Thiel et al. 
8| have studied the transport properties of the LaAlOs/SrTiOa interface using the electric 
field-effect. Measurements on samples without indications of oxygen defects show that the 
LaAlOs surface layer must reach a critical thickness of 4 unit cells for the interface to be 
conducting. 

From the theoretical point of view, charge compensation at LaAlOs/SrTiOs interfaces 
has been investigated by Park et al. Qj. Electron doping of Ti atoms in combination with 
the Jahn- Teller effect leads to metallicity in their band structure calculations, based on the 
local density approximation (LDA). In addition, applying the LDA+U scheme, Pentcheva 
and Pickett [ICl] have found ferromagnetically aligned spins due to occupied d^y orbitals in a 



checkerboard of Ti'^"*" sites. Beyond structural relaxation effects at the 



ous modifications of the electronic structure are caused by surfaces H 



leterointerface, seri- 



m. FortheLaAlOa 



(001) surface, atomic reconfiguration due to the polarity discontinuity at the vacuum inter- 
face is accompanied by a remarkable redistribution of the surface electron density, see Lanier 
et al. [l3| and the references therein. In spite of the relevance of these results for heteroint- 
erfaces near surfaces, first principles band structure calculations accounting for the vacuum 
interface are missing so far. Our present work therefore addresses LaAlOs/SrTiOa multi- 
layer configurations terminated by an (A102)~ surface. We discuss the interplay between the 
electronic states at the heterointerface and their distance from the surface. Modifying this 



distance turns out to be a most useful instrument for systematically tuning the transport 
properties. 

The band structure results discussed in the following rely on the scalar-relativistic aug- 
mented spherical wave (ASW) approach [ij]. We apply a recently implemented code which 
accounts for the non-spherical contributions to the charge density inside the atomic spheres 
llSl . The ASW method is advantageous for dealing with complex structures comprising a 
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large number of atoms [16|, 117|, [181]. We model heterointerfaces by long tetragonal supercells 
based on the parent perovskite structure. The principal axis of these cells gives rise the crys- 
tallographical c-axis, which runs perpendicular to the surface. To be specific, a supercell 
comprises a layer of 4 SrTiOs unit cells sandwiched between layers of 3-4 LaAlOs unit cells, 
compare figure [H The LaAlOa layers themselves are terminated by a vacuum layer of at 
least 12 A thickness, allowing us to apply periodic boundary conditions. As superstructures 
perpendicular to the interface have not been reported, they are not taken into account. 
Moreover, we use the symmetric crystal structure of high temperature cubic SrTiOs for the 
whole supercell. Lattice mismatch between SrTiOs and LaAlOs layers has to be avoided via 
a common lattice constant. Since LaAlOa is grown on a SrTiOs substrate in the experiment, 
it is reasonable to choose the measured value for cubic SrTiOs, which amounts to 3.91 A 
[l9(]. In the LaAlOs region, this approximation artificially elongates the lattice constant by 



(less than) 2.5% with respect to the experimental value of 3.81 A 19|]. Our lattice setup is 
supported by experimental data as well as theoretical structure optimization. 

Structural relaxation at both the LaAlOs/SrTiOs interface and the surface is taken into 
account via a minimization of the atomic forces. For this purpose we use the Wien2k program 
package [20|. This is a famous full-potential hnearized augmented plane wave code having 



shown great capability in the optimization of surfaces as well as interfaces 2l|, |22|. The 
surface relaxation starts out from a supercell consisting of two LaAlOs unit cells along the 
c-axis, terminated by a vacuum layer. As compared to the bulk atomic positions, the surface 
(A102)~ layer developes a distinct buckling, because the O atoms are shifted to the vacuum 
by 0.02 A and the Al atoms move 0.05 A in the opposite direction. In the adjacent (LaO)^ 
layer, the La atoms follow the surface O sites by 0.03 A, while the O atoms retain their 
bulk positions. Turning to the relaxation effects at the heterointerface, we use a supercell 
of touching LaAlOs and SrTiOs layers, each covering 4 perovskite unit cells in c-direction. 
In the (LaO)"*" layer at the interface, only the O atoms move 0.07 A off the contact, pushing 
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TABLE L Metal-0 bond lengths at the surface and the interface. The names of the O sites refer 
to the layer they are located in. For the Sr-0 bond length the short/long bond points off/towards 
the interface. The structure optimization has started from bond lengths of 1.95 A (Ti/Al-0) and 
2.76 A (Sr/La-0), while the intraplanar metal-0 distances keep their original values. 

neighbouring Al atoms 0.02 A in the same direction. In contrast, structural modifications 
are more pronounced in the (Ti02)° interface layer: Whereas the Ti atoms are shifted off 
the contact by 0.04 A, the O atoms approach it by 0.06 A. Finally, the neighbouring (SrO)° 
layer reveals shifts of the Sr and O atoms by 0.03 A off and towards the contact, respectively. 
These findings as well as the optimized bond lengths summarized in Table |T] agree well with 
the structural data reported by Park et al. l|. 

For the structure optimization we use the generalized gradient approximation with a 
mixed linear-augmented-plane-wave and augmented-plane-wave plus local-orbitals basis. In 
the surface and contact calculation the charge density is represented via about 4,000 and 
11,000 plane waves, respectively, and the k-space grid has 10 and 15 k-points in the irre- 
ducible wedge of the Brillouin zone. In addition, the Perdew-Burke-Ernzernhof parametriza- 
tion of the exchange-correlation potential is applied. On the contrary, the ASW supercell 
calculation uses the local density approximation within the Vosko-Wilk-Nusair scheme. Aim- 
ing at a correct representation of the crystal potential in large voids of the supercell and 
near the surface, the physical atomic spheres are complemented by additional augmentation 
spheres at carefully selected interstitial sites. In the case that the surface layer comprises 
3 and 4 LaAlOs unit cells, respectively, it is sufficient to dispose 153 and 173 additional 
spheres in order to keep the linear overlap of the physical spheres below 17% and of any 
pair of spheres below 22%. Along with 58 and 68 physical spheres, the supercells entering 



(a) Thickness of the LaAlOa surface layer: 3 unit cells 
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(b) Thickness of the LaAlOs surface layer: 4 unit cells 
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FIG. 3: (Color online) Band structure and partial Ti 3d and O 2p DOS (per atomic site) for the 
interface (Ti02)'' layer, which is separated by (a) 3 and (b) 4 LaAlOs unit cells from the surface, 
respectively. The weighted bands refer to the simple cubic Brillouin zone of the parent perovskite 
structure, showing only states with Ti 3d contributions beyond a certain threshold. 

the ASW calculation thus comprise 211 and 241 augmentation spheres in total. The basis 
set taken into account in the secular matrix in each case consists of La 6s, 6p, 5d, 4/, Sr 5s, 
5p, 4d, Ti 4s, 4p, 3d, and O 2s, 2p orbitals, as well as states of the additional augmentation 
spheres. During the course of the band calculation the Brillouin zone is sampled with an 
increasing number of up to 110 k-points in the irreducible wedge, which ensures convergence 
of the findings with respect to the fineness of the k-space grid. 

To start the discussion of the electronic structure results, we first analyze bulk SrTiOs for 
comparison. The band structure is depicted in figure [2] along selected high symmetry lines 
in the first Brillouin zone of the simple cubic lattice, where the symmetry points are defined 



by standard reciprocal vectors: T = (0, 0, 0), M = (|, i, 0), X = (0, |, 0), and Z = (0, 0, |). 
Contributions of Ti 3d t2g states are represented by the length of the bars added to every 
band and k-point. Since the system is insulating, the conduction band is empty, with the 
minimum at the F-point. One of the bands giving rise to the conduction band minimum has 
hardly any dispersion along the lines T-X and T-Z, which are degenerate by symmetry, thus 
along the cubic lattice vectors. The remaining states are characterized by strong dispersion, 
reflected by the total band width of about 2.3 eV. Corresponding densities of states (DOS) 
are given on the right hand side of figure [2l In the energy range shown, contributions of Sr 
states almost vanish. As to be expected from a simple molecular orbital picture, the occupied 
and unoccupied group of bands is dominated by O 2p and Ti 3d contributions, respectively, 
since orbital overlap leads to bonding and antibonding molecular states. Admixtures in the 
energy range dominated by the respective other states are due to TiO-hybridization, and 
thus are observed mainly for cr-type overlap. These findings agree well with previous band 
structure data, see (23[, for example. 

Turning to the LaAlOs/SrTiOs interface, band structures and DOS curves are given in 
figure [3] in the same representation as used for bulk SrTiOs, see figure El In particular, 
the weighted bands refer to the simple cubic lattice of the parent perovskite structure. For 
clarity, we show only such states with sufficient Ti 3d contributions from the electron-doped 
(Ti02)° layer right at the interface. Figure [3] compares the data obtained for different 
surface-interface distances: The thickness of the LaAlOs surface layer amounts to (a) 3 and 
(b) 4 perovskite unit cells. As concerns the gross features of the band structure as well as the 
DOS, the contact Ti 3d states largely resemble the bulk SrTiOs findings. However, the group 
of conduction bands now is filled partially, which leaves the interface in a conducting state. 
To be more specific, dispersion disappears completely along the line F-Z due to the spacial 
restriction induced by the surface. Moreover, since the single occupied Ti 3d band shows 
strong dispersion within the a6-plane, it belongs to nearly ideal two-dimensional states. 
They are visible in the DOS just below the Fermi energy. 

Our results confirm the observation of Thiel et al. J8| that a growing thickness of the 
LaAlOs surface layer enhances the conductivity of the heterointerface. Because surface elec- 
tronic states leak out into the vacuum area, the charge carrier density is reduced in the 
LaAlOs surface layer. As a consequence, less charge is available for doping the interface. 
When the distance between interface and surface increases, the charge depletion declines 



because of electronic screening. This trend is obvious in figure [31 For the thinner LaAlOa 
surface layer the two-dimensional band reaches down to — 0.5eV, whereas the minimum is 
found at — 0.7eV in the case of the thicker layer. The occupation therefore increases system- 
atically with the number of LaAlOs cells separating the vacuum from the hetero interface. 

As concerns the qualitative dependence on the thickness of the LaAlOs layer, theory and 
experiment agree excellent, while a critical number of 4 unit cells for conductivity is not 
reproduced by the calculation. Instead, conductivity is likewise achieved in the case of 3 
unit cells, and even for thinner surface layers, as checked by additional calculations. This 
discrepancy probably is connected to the underestimation of the band gap in pure SrTiOa, 
tracing back to the local density approximation. In figure [H the gap amounts to some 
2.1 eV, which is significantly less than the experimental value 3.2 eV. Similar shortcomings 



are known for Ti02 2J] and LaAlOs [25]. As a consequence of the underestimated band 
gap, the occupation of the two-dimensional band at the heterointerface is overestimated. 
The transition from insulating to conducting behaviour therefore cannot not be caught 
quantitatively. However, these shortcomings do not affect our previous conclusions about 
surface effects on heterointerfaces. In fact, the dependence of the interface charging on 
the surface layer thickness is hardly affected by the structural relaxation, since largely the 
same behaviour is found for non-relaxed interfaces. This strongly supports an interpretation 
in terms of charge depletion due to the proximity to the surface. In addition, electronic 
correlations would not be expected to be critical for the mechanism under consideration. 
Our results thus do not depend on the specific compounds giving rise to the interface. 

In conclusion, we have presented first-principles band structure results for LaAlOs/SrTiOs 
heterointerfaces in the vicinity of a surface. As a consequence of electron-doping due to 
charge transfer from the LaAlOa surface layer, we find the interface (Ti02)° layer to be 
conducting. As the dispersion of the conduction band is strong and reveals an ideal two- 
dimensional character, a quasi two-dimensional gas of very mobile electrons is formed. In the 
case of an ultrathin LaAlOs layer, surface effects are of special importance for the electronic 
structure at the heterointerface. Since electronic states leak out into the vacuum, charge 
is distributed close to a surface. For the heterointerface, this implies a rigid shift of the 
electronic bands because of a modified Fermi level. In other words, charge transfer towards 
a surface counteracts an intrinsic interface doping. Since this effect declines with the surface- 
interface distance, the charge carrier density can be adjusted via the thickness of the surface 

8 



layer. 
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